Classical and quantum correlations under decoherence 
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Recently some authors have pointed out that there exist nonclassical correlations which are more 
general, and possibly more fundamental, than entanglement. For these general quantum correlations 
and their classical counterparts, under the action of decoherence, we identify three general types of 
dynamics that include a peculiar sudden change in their decay rates. We show that, under suitable 
conditions, the classical correlation is unaffected by decoherence. Such dynamic behavior suggests 
an operational measure of both classical and quantum correlations that can be computed without 
any extremization procedure. 

PACS numbers: 03.65.Ta, 03.67.-a, 03.65.Yz 



It is largely accepted that quantum mutual information 
is the information-theoretic measure of the total correla- 
tion in a bipartite quantum state. Groisman et al. 
inspired by Landauer's erasure principle 0, gave an op- 
erational definition of correlations based on the amount 
of noise required to destroy them. From this definition, 
they proved that the total amount of correlation in any 
bipartite quantum state (pab) is equal to the quantum 
mutual information [I(pa-.b) — S(pa) + &(Pb) ~ S(pab), 
where S(p) = — Tr(plog 2 p) is the von Neumann entropy 
and Pa(b) — ^b{A){pab) is the reduced density oper- 
ator of the partition A(B)\. Another argument in fa- 
vor of the claim that quantum mutual information is a 
measure of the total correlation in a bipartite quantum 
state was given by Schumacher and Westmoreland 
They showed that, if Alice and Bob share a correlated 
composite quantum system that is used as the key for 
a "one-time pad cryptographic system" , the maximum 
amount of information that Alice can send securely to 
Bob is the quantum mutual information of the shared 
correlated state. 

We are interested here in the dynamics of both quan- 
tum and classical correlations under the action of noisy 
environments. For these purposes, it is reasonable to as- 
sume that the total correlation contained in a bipartite 
quantum state may be separated as J(pa-.b) — Q{pab) + 
C{pab), owing to the distinct nature of quantum (Q) and 
classical (C) correlations [l], 0, 0, 0] ■ Some proposals for 
characterization and quantification of Q and C in a com- 
posite quantum state have appeared in the last few years 
jE 0, IE 0, 0, 0j- The quantum correlation, Q(pab), 
between partitions A and B of a composite state can 
be quantified by the so-called quantum discord, V{pab), 
introduced by Ollivier and Zurek Q. Such a quantum 
correlation is more general than entanglement, in the 
sense that separable mixed states can have a nonclas- 
sical correlation that leads to a nonzero discord. It mea- 
sures general nonclassical correlations, including entan- 



glement. For separable mixed states (unentangled states) 
with nonzero discord, this quantum correlation provides 
a speed up, in performing some tasks, over the best 
known classical counterpart, as was shown theoretically 
f]~oj ] and experimentally [ill ] in a non-universal model of 
quantum computation. Therefore, such a nonclassical 
correlation might have a significant role in quantum in- 
formation protocols. For pure states, we have a special 
situation where the quantum correlation is equal to the 
entropy of entanglement and also to the classical corre- 
lation. In other words, Q(pab) — C(pab) — Z(pa-.b)/2 
[3, 0]. In this case, the total amount of quantum cor- 
relation is captured by an entanglement measure. On 
the other hand, for mixed states, the entanglement is 
only a part of this more general nonclassical correlation, 
Q(pab) 0,113, 1! I- A quantum composite state may also 
have a classical correlation, C(pab), which for bipartite 
quantum states can be quantified via the measure pro- 
posed by Henderson and one of us Since we assume 
that the total correlation is given by the quantum mu- 
tual information and if we adopt the definition of classical 
correlation given in Q(pab) turns out to be identi- 
cal to the definition of quantum discord in Ref. 0]; in 
other words, Q{pab) = V(pab) = Z(pa-.b) - C(pab), as 
already noted in Ref. |I3 |. 

We have identified three different kinds of dynamic 
behavior of C and Q under decoherence, which depend 
on the "geometry" of the initial composite state and on 
the noise channel: (i) C remains constant and Q de- 
cays monotonically over time; (it) C suffers a sudden 
change in behavior, decaying monotonically until a spe- 
cific parametrized time, psc (to be defined below), and 
remaining constant thereafter, while Q has an abrupt 
change in its rate of decay at psc, becoming greater than 
C within certain parametrized time interval; and (iii) C 
and Q decay monotonically. For two-qubit states with 
maximally mixed marginals we show which conditions 
lead to the different types of dynamic behavior, for cer- 
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tain noise channels (i.e., phase flip, bit flip, and bit-phase 
flip). We also recognize a symmetry among these chan- 
nels and provide a necessary condition for C to remain 
constant under decoherence, which enables us to define 
an operational measure for both classical and quantum 
correlations. 

Let us start with the definition of classical correlation 



4] = 



C(pab) = max [S{p A ) ~ S {Uj }{pA\B)] 



(i) 



no- 



where the maximum is taken over the set of 
jective measurements {II? } [l3| on subsystem B [3], 

S{iij}(P a\b) = {j^a) i s t ne conditional entropy 

of subsystem A, given the knowledge (measure) of the 
state of subsystem B, p A = Tr B (JljPAB^j)/ qj, and 
qj = Tr AB (p AB Tlj). 

We consider the scenario of two qubits under local de- 
coherence channels. The evolved state of such a sys- 
tem under local environments may be described as a 
completely positive trace preserving map, e(-), which, 
written in the operator-sum representation, is given by 

MM 



£ ( MB ) = £r^rfW 



,(S)f r (A)t 



where Tr°' (k = A, B) are the Kraus operators that de- 
scribe the noise channels A and B. 

For simplicity, let us consider a class of states with 
maximally mixed marginals (pa(B) = 1a(b)/2), de- 
scribed by 



PAB 




(2) 



where a\ is the standard Pauli operator in direction i 
acting on the subspace k = A,B, ci € 3? such that < 
\ci\ < 1 for i = l,2, 3, and 1a(b) is the identity operator 
in subspace A(B). The state in Eq. ((2]) represents a 
considerable class of states including the Werner (|ci| = 
|c 2 ] = |c 3 | = c) and Bell (|ci| = |c 2 | = |c 3 | = 1) basis 
states. 

Phase flip channel. This is a quantum noise pro- 
cess with loss of quantum information without loss of 
energy. For this channel, the Kraus operators are 
given by r< A) = diag(y/l - p A /2, y/l- PA /2) ® 

Is, r[ A) = diag{^/p~^,-^/p~I/2) ® l s , T[ B) = 
1 A ® d iagj^ /l - pbT^ , y/l-p B /2), and t[ B) = 1 A <g> 
diag(y/pB~/2, — \/pb/2), written in the subsystem basis 
{|0}fc,|l)fe}, k = A,B. We are using p A(B) (0 < p A {B) < 
1) as parametrized time in channel A(B). We consider 
here the symmetric situation in which the decoherence 
rate is equal in both channels, so p A = ps = P- 



The description of the dynamical evolution of the sys- 
tem under the action of a decoherence channel using the 
parametrized time p is more general than that using a 
specific functional dependence on time t, in the sense 
that it accounts for a large range of physical scenarios. 
For example, for the phase damping channel (the phase 
damping and phase flip channels are the same quantum 
operation (Til]), we have p = 1 — cxp(— jt), where 7 is the 
phase damping rate fl7j |. 

The density operator in Eq. ([2]) under the multimodc 
noise channel, s(p AB ), has the eigenvalue spectrum: 



Ai = ~ [1 - a - p - 7] 
A3 = \ [1 + a - p + 7] 



A 2 = - [1 - a + + 7] , 

A 4 = ~[l + a + /3-7]> (3) 



with a = (1 — p) 2 Ci, /3 = (1 — p) 2 C2, 7 = C3, and the von 
Neumann entropies of the marginal states remain con- 
stant under phase flip for any p, S \Tia(b) £ (pab)] = !• 
To compute the classical correlation (TTJ) under phase 
flip, we take the complete set of orthonormal projec- 
tors {IL; = |6,) (%| ,j =||, -L}, where |6||) = cos(0) |0) + 



e^sin((9) |1) and |6j 



^sin(0) |O)-cos(0) |1). Then 



the reduced measured density operator of subsystem A 
under phase flip, J? A = Tr# \Tljs(p A B)Tlj]/ qj, will have 
the following eigenvalue spectrum: 



^ = i{2±[2 7 2 + a 



(3 2 + (2 7 2 -a 2 - f3 2 ) cos (40) 
+2{a 2 - /3 2 ) cos (20) sin 2 (20)] 1/2 }, (4) 
and qj = 1/2, for j =||, _L. From Eq. ([T]), it follows that 



C[e(p AB )] = l-mm[s (p A ) 



(5) 



since = £^ and hence S (pI) — S (p A )- The classi- 
cal correlation and the quantum correlation under phase 
flip may be written, respectively, as 

C [e (pab)] = E 1 + ( 2 ^ 1 °S 2 ( 1 + M 
fc=i 

4 

Q [e {p AB )] = 2 + Afc log 2 A fe - C [e(p AB )] , (6b) 
fe=i 

where \ = max (|a| , , I7I), which depends on the re- 
lation between the coefficients Cj in state <j2j) and on the 
parametrized time p. 

(i) If l c 3| > l c i| > l c 2| m (J2), the minimum in ([5]) is ob- 
tained by = 4> = 0. The classical and the quantum cor- 
relations under phase flip will be given in Eqs. (|6a[) and 
(|6b|) . respectively, with \ = \ c z\- I n this case, the classi- 
cal correlation C [e(p A B)\ is constant (it does not depend 
on the parametrized time p) and equal to the mutual 
information of the completely decohered state (p = 1), 
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C{pab) = C [e(pab)\ 
turn correlation [Eq. 



= 1 

Mb 



£{PA-.b)\ 



while the quan- 



decays monotonically. 
(ii) If |ci| > |c 2 |,|c 3 | or |c 2 | > |ci|,|c 3 |; and |c 3 ^ 
0, we have a peculiar dynamics with a sudden change 
in behavior. C decays monotonically until a specific 
parametrized time, psc = 1 — y |c3|/max(|ci| , \c2\), 
and from then on C remains constant. For p < psc, 
the minimum in (O is achieved when 9 = 7r/4, — 
(if |ci| > |c 2 |) or = ir/2 (if |d| < |c 2 |), and X = 
(1 — p) 2 max(|ci| , |c 2 |). Thus, C decays monotonically. 
On the other hand, for p > psc, the choice 9 = 4> = 
leads to the minimum in ((5]) and x = l c 3|- Then C sud- 



denly becomes constant at p = pscs C 



£(PA:B)\ p=1 



~{pab) 



p>psc 



and the decay rate of Q changes sud- 
denly at p = psc- In Fig- 1> we depict this peculiar 
behavior for a given choice of parameters and, in Fig. 2, 
we show the values of the sudden change parametrized 
time, p sc , as a function of Ci and c 2 . 

(iii) Finally, if | C3 1 = 0, we have a monotonic decay of 
both correlations C and Q. 



that the quantum correlation may be greater than the 
classical one for some states, for example e(pa.b)\ v — vso - 
It is worth mentioning that this peculiar sudden change 
in behavior is a different phenomenon from entanglement 
sudden death Indeed, it seems that these 



correlations do not present sudden death [21 1 . 





FIG. 1: Classical C [e(pab)] (dashed line), quantum 
Q[e(Pab)\ (solid line), and total I [e(pa-.b)] (dotted line) 
correlations under phase flip. We have set, in this figure, 
ci = 0.06, c 2 = 0.42, and c 3 = 0.30 . For this state the 
sudden change occurs at p sc — 0.15, and Q is greater than 
C for 0.09 < p < 0.20. At p = 0.09 and p = 0.20, we have 
Q [e(pab)] =C[e (pab)] = I [e(pa-.b)]/ 2, as happens for pure 
states. 



The dynamic behavior of correlations under the phase 
flip channel described in Fig. 1 is quite general. Such a 
sudden change in behavior occurs also when we consider 
the bit flip and the bit-phase flip channels [of course un- 
der other conditions on the Ck 's in state J^J]. Moreover, 
these results contradict the early conjecture that C> Q 
for any quantum state [j], 0, EH- Here, we have shown 



FIG. 2: Sudden change parametrized time. p 3C as a function 
of c\ and C2, for C3 = 0.1, under a phase flip channel. In the 
regions where p sc = or p ac = 1 there is no sudden change. 



Bit flip channel. The Kraus operators are [15 



<A) 



p/2a { x 



and T 



(B) 



diag{y r T- 

(B) 
LB, X = 



p/2,yr^/2) <g> 1 B , ri- 



(A) 



T"" 1 = l A ®diag{ v r l 



p/2, s/T^p/2), 
The eigenvalue spectrum of 
e (pab) is given by where the variables now take the 



form a = c±, /3 = (1 — p) c 2 , and 7 = (1 — p) C3. The 
correlations can again be written as ([oa]) and (j6b[) . The 
dynamic behavior of C and Q under bit flip is symmetrical 
to that for the phase flip channel (just exchanging c\ and 
C3). Type {%) dynamics is obtained when |ci| > |c 2 | , | C3 1 . 
Type (ii) occurs for |c 3 | > |ci|,|c 2 | or |c 2 | > |ci|,|c 3 |, 
and |ci| 0, with a sudden change in behavior of C and 
Q at psc = 1 - \/\ci\/ max(|c 2 | , |c 3 |). Finally, if \a\ = 0, 
we have type (iii) dynamics. 

Bit-phase flip channel. Now, the Kraus operators are 
T { Q A) = diag(y^pj2, y^p~/2) ® 1 B , r[ A) = 




andT 



T ( B) = l A ®diag(y / l-p/2, y/l - p/2), 



'1b 

, - lA®^/pj2ay B> . The variables in Eq. ([3]) turn 
out to be a = (1 —p) c±, j3 = c 2 , and 7= (1 — p) c 3 . C 
and Q under bit-phase flip can again be written as |6a[) 
and (j6bj) . respectively. Once more, the conditions for the 
various types of dynamics are obtained by swapping c 2 
and c 3 in the phase flip channel. For type (ii) dynamics, 
we now have psc = 1 - y 1 |c 2 |/ max(|ci | , \c 3 \). 

Necessary conditions for C to remain constant under 
decoherence are the following: 



n r (s) 

L1 J ' 1 k 



= 0, V j, k. 



(7) 
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These relations depend on the angles 9 and <f> that define 
the minimum in (HI). For the channels mentioned above, 



r 



(B, 



1 B and T\ 



(B) 



with i = 1 for the bit flip, 



i = 2 for the bit-phase flip, and i = 3 for the phase flip. 
Hence, condition (0 will be satisfied when the projective 
measurements that reach the minimum in Eq. ([5]), IP,-, 

( B) | A. 

are performed on eigenstates of <j\ 22]. On the other 
hand, the angles 6 and <f> that define the minimum in 
Eq. (O depend on the "geometry" of the initial state. 
When the larger component of state in Eq. ((2|) is in the 
direction 1, 2, or 3, C remains constant under bit flip, 
bit-phase flip or phase flip, respectively. 

The fact that, for a given state, the classical correlation 
can remain unaffected by a suitable choice of noise chan- 
nel, e, immediately suggests an operational way (with- 
out any extremization procedure) of computing classical 
and quantum correlations. It could be done as follows: 
depending on the state "geometry" , we send its compo- 
nent parts through local channels that preserve its clas- 
sical correlation, so that the quantum correlation will be 
given simply by the difference between the state mutual 
information X(pa.b) and the completely decohered mu- 



tual information. X 



£{pa-.b)\ p=1 



Q(pab) = X(pa-.b) - X £{pa-.b)\ p=1 
since X(p A . B ) = Q(pab) + C(pab) and 
C(pab) =X e(pA:s)L = 



A suitable channel for the class of states described by Eq. 
([2|) is chosen which satisfies condition ((7|) as discussed 
above. 

A problem to be addressed before such a measure can 
be used for a general state is to establish a protocol to 
find the map (if this map exists) which leaves the classical 
correlation unaffected [23j. This suggests an interesting 
research program to develop an operational way of inves- 
tigating the role of quantum and classical correlations in 
many scenarios, such as quantum phase transitions 241 ] . 
non-equilibrium thermodynamics [251 ]. etc. 
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